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ABSTRACT 


This report is a detailed instruction manual for users 
of the 1962 LASL lens design code, which is a general-purpose 
computing program for the evaluation and automatic design of 
optical lens and/or mirror systems with surfaces generated by 
conic sections. Described in relation to the program are IBM- 
7090 machine operations, prescription parameters, ray-tracing 
parameters, weights on lens performance, lens-parameter sub- 
stitutions, increments to lens parameters, designing instruc- 
tions, and parameter interactions. The sample calculation 
for a Lister-type lens is given in great detail. The general 
procedure used to design a special-purpose zoom lens is given, 
along with the final lens prescription and the performance 


characteristics. 
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INTRODUCTION 


The purpose of this paper is to explain the use of the 
1962 LASL automatic lens-designing code # and so to make optical 
lens designing available to any engineer or scientist with a 
modest knowledge of optics. The code gives to problems of lens 
design a general automated solution which makes full use of the 
digital computer's capabilities and abandons the classical 
approach requiring direction by the skilled specialist at every 
stage. 

To study and design lenses, the code makes use of geomet- 
rical optics, the geometrical analysis of the light rays passing 
through the lens. Three laws of geometrical optics, those 
governing rectilinear propagation, refraction, and reflection, 
are used. By using these laws to trace light rays through a 
lens, we can analyze lens performance. The main criterion for 
lens performance is the image-spot size at the focal plane of 
interest. While image defects are analyzed precisely by eval- 
uating the sizes and positions of image spots at the focal sur- 
face, the traditional Seidel aberrations are not evaluated. In 


specifying, selecting, or designing lenses for particular tasks, 


this ability to analyze lens performance is invaluable. 

The 1962 code developed by C. A. Lehman for the IBM-7090 
Data Processing System is an adaptation of the automatic, 
iterative program developed in 1959 by J. C. Holladay for the 
TBM To.© The LASL history of lens designing with computers 
began in 1950 when R. H. Stark applied floating-decimal 
computer notation to the LASL IBM-CPC computer .2 Then 
W. A. Allen and Stark’ made the first application to lens 
designing when they produced an easy to use ray-tracing program. 
In the same year R. E. von Holdt improved this program by 
substituting vectorial technique? for the algebraic methods 
used by Allen and Stark. These programs, however, limited 
themselves to ray tracing and required interpretation by an 
expert. Until 1959, attempts to broaden the system on the CPC 
and the Maniac met with indifferent success. 

The Holladay program made statistical analyses of lens 
performance from skew ray traces. It determined differences 
in lens performance resulting from design-parameter alterations. 
Using these differences, it computed the linear combination 
which would give an improved lens prescription. With this 
system it was possible to design a complete lens by means of 
the computer. The success of this program created a desire 
for a code which would be both easier to use and able to 


solve more complex problems, 


Consequently, Lehman has developed this faster code for 
the IBM-7090 Data Processing System. This new code can produce 
a lens design which is the result of a simultaneous analysis 
of as many as five different conditions of use. (A specific 
example is the zoom lens--to be described later--with five 
different focal lengths, apertures, and angles of view.) This 
code, in addition, can trace rays through cylindrical, displaced, 
and tilted surfaces, a feat not previously possible. Also new 
is simultaneous design for a maximum of six wave lengths. Use 
of the code is facilitated by the following automation improve- 
ments: design-parameter cycling, determinant adjustment, 
negative thickness prevention, and ray vignette informing. 

The theory of this lens-designing program is treated in 
detail in Lehman's report and will only be summarized in the 
following two sections. All information about lens performance 
is derived from rays traced by a skew ray-tracing routine. 

The deviation of the traced rays from their specified image 
points is an error that is reduced by alteration of the lens 
parameters as directed by the code's designing routine, which 
uses the least-squares technique. Since the subroutines for 
obtaining the square root and the linear matrix solution are 
well-known mathematical operations, they are not described. 


Also not described are the subroutines for on-line card reading, 


on-line and off-line printing, punching binary prescription- 
data cards, and punching decimal point-plot cards. 


RAY TRACING 


Skew rays are traced through the lens or mirror system 
by the methods of geometrical optics. Computation is performed 
by the vectorial technique. The ray is traced in two stages 
to avoid the loss of accuracy which may occur in the machine 
solution of a quadratic equation when there is a long optical 
path. First a linear equation is used to locate the ray on the 
vertex tangent plane of the optical surface. Next the ray is 
traced to the optical surface and refracted or reflected. Then 
the unit vector for the refracted or reflected ray is calculated. 
This cycle of operations is repeated for each optical surface. 
Finally, the intersection of the ray with the focal surface (or 


other aiming point) is computed. 


DESIGNING 


Each ray that is traced through the lens system would 
intersect the specified image point if the lens gave perfect 
performance. In an imperfect lens the deviation of the ray 
from the specified image point is an error which can be weighed 


to make a weighed-error component. The sum of the weighed-error 
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components of all the traced rays is an error vector which is 
a measure of lens performance. To optimize the lens performance 
it is necessary to minimize the magnitude of the error vector 
by alteration of the lens parameters. The parameter alterations 
necessary to effect an improvement are determined by a least- 
squares method which is here applied by an arbitrary approxi- 
mation of Newton's method of solving for minima of complicated 
functions. Selected design parameters of the lens are individ- 
ually altered by a small amount, and the error vectors are 
computed for each. These error vectors are used to generate 
a matrix which is solved linearly to obtain an approximation 
of the parameter changes necessary to reduce the error vector 
and hence improve the performance of the lens. 

Since the above solution is only an approximation of the 
desired error-vector minimum, only a modest parameter change 
is allowable, a maximum of perhaps 10 to 100 times the altera- 
tions used to make the computations. A maximum of 10 parameter 
increments can be calculated simultaneously, but experience 
indicates that on the preliminary design fewer should be used 
because of the increased dependence of the variables upon each 
other as the number increases. An iterative procedure is used 
to arrive at the minimum error vector obtainable with a specific 


set of parameter alterations. 
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WHAT THE PROGRAM CAN DO 


The 1962 LASL lens designing code is a general-purpose 
computing program for the investigation of lens systems with 
as many as 98 surfaces. The program can determine the per- 
formance of an existing lens design. It can automatically 
improve an existing lens design to meet the conditions required 
by specific problems. And, finally, it can automatically obtain 
the optimum lens design from a start consisting of any desired 
arrangement of refracting and/or reflecting surfaces. 

The lens or mirror surfaces most used are either spheres 
or cylinders but surfaces with any conic section may be used. 
The cylinder axis is always parallel to the x axis. Each of 
the surfaces of revolution may be displaced or tilted as desired. 
There is a wide choice possible in the number and distribution 
of the rays which may be traced through the lens. The positions 
and sizes of the entrance and exit pupils are readily controlled. 
The lens focus may be fitted to a specified focal surface, or the 
shape of the focal surface may be determined. The machine can 
design simultaneously for as many as seven object points, six 
wave lengths, and five different conditions of use. 

The description which follows is a graphically illustrated 
exposition aimed at aiding the engineer to master the complex- 


ities of the use of the code. It has many numerical examples 


Ue 


of the data and of the program instructions which have been 


found useful for designing. 
MACHINE OPERATION 


This code is for an IBM-7090 Data Processing Gyacene with 
the following components: a 7302 Core Storage (32,768-word), 

a 7606 Multiplexor, a 7607 Data Channel, SHARE II control 
panels, a 711-II Card Reader, a 716 Printer, a 721 Card Punch, 
a 729-IV Magnetic Tape Unit connected as logical tape A3, and 
a 1401 Data Processing System to interpret the tape record. 

All data are read from punched cards and must be addressed to 
cells specifically reserved for each item so that the operation 
programs can locate the specific data when needed. 

The first operation which the operator must perform is to 
load the 213-card code into the card reader, where the program's 
coded instructions are read and sent to the core memory. This 
operation is followed by loading the lens data, which are on 
binary and/or decimal punched cards (columns 1 to 72). The lens 
data include the lens prescription and the parameters for ray 
tracing and designing. There are 10 individual program-control 
instructions on decimal punched cards, which are used to select 
the sequence of machine calculations. 


All cards with decimal data and all operation instructions 
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are printed on line as soon as read. The tabulated lens data 
will be printed on line if the program discovers an error in 
the data. After completion of the calculations ordered by 
each program-control instruction, the specified answers and 
the latest lens prescription are recorded on the A3 tape. An 
on-line print of the tape record will also be obtained when 
sense switch six is down. The 721 Card Punch will make binary 
punched cards with the latest complete lens data (instruction 
*6 below [plus *3 and *4]) if sense switch one is down at the 
end of any of the ten routine program calculations. These lens 
data on binary punched cards are used to restart the problem 
for the next machine run. Lens-data changes can be read in 
before the start of any of the routine program calculations. 
The 1401 Processing Unit will, if ordered, supply the spot- 
diagram coordinates on decimal punched cards at the same time 
that the off-line print is made. 

Decimal numerical (and alphabetic) data are stored in the 
machine's magnetic core memory after being read from punched 
cards and translated to binary notation. The cards are punched 
according to the conventions outlined below. 

i. Integers prefixed by L or L- order the program to store 
the next number at that address (L-200 to L1999). Succeeding 
numbers go to succeeding addresses. 


2. Integers prefixed by C order the program to clear that 
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many addresses in sequence. 

3. Integers prefixed by X or X- are tally numbers. 

4. Integers prefixed by * are program-control instructions. 

5. Decimal numbers prefixed by + or - are positive or negative 
numbers representing the parameters. These numbers may contain 

as many as 10 digits and a decimal point. They may each be 

followed by a power-of-ten multiplier consisting of an integer 

(1 to 38) prefixed by E or E-. 

6. The problem identification (any combination of machine 

characters and spaces heading each prescription print) placed 

between $'s will be stored as written (six per address) if 

prefixed by H and addressed to L-200 through L-171. The char- * 500 
acters in L-200 and L-199 also appear on the print and the 7a 
punched cards for spot-diagram plotting. 

The program reads decimal punched cards automatically. 
Specified answers to calculations are printed automatically. 

The program control instructions and the calculations ordered 
by each are described below. 

The symbol *1 instructs the program to read data on binary 
punched cards until a blank card is encountered, at which time 
the program returns to reading decimal punched cards. 

*2 instructs the program to proceed with the lens designing 
(cycles and specifications scheduled at L-79 through L-65, L-49 


+ 
Marginal numbers are an index of storage addresses. 
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through L-3, and L1010 through L1499) and tabulate answers 


showing what has been accomplished. 


*3 instructs the program to determine lens performance by 


multiple ray tracing (scheduled at L-170 through TOS a= 59 
through L-56, and L1500 through L1999) and tabulate the answers 
giving the sizes and positions of the images. 

*4 instructs the program to trace the upper and lower 


rim rays from the maximum object height and to tabulate the 


ray position and direction vectors and the distance from the 


preceding surface at each optical surface. 


*5 instructs the program to trace paraxial and rim rays 


from zero object height and to tabulate the ray position and 
direction vectors and the distance from the preceding surface 


at each optical surface. 


*6 instructs the program to make binary punched cards with 


the latest lens data. 


*7 instructs the program to calculate and tabulate the 


scaled ray-intersection coordinates for spot-diagram plotting. 


*8 instructs the program to calculate, tabulate, and make 
decimal punched cards with the scaled ray-intersection coordi- 
nates for spot-diagram plotting. 

*Q instructs the program to tabulate the latest lens 
prescription. 


*10 instructs the program to tabulate all of the lens data 
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that are in storage. 


*O instructs the program to stop, signifying the end of 


the run. 


PRESCRIPTION PARAMETERS 


To determine the performance of an existing lens pre- 
scription we must know or find the wave-length (color) range, 
the distance to the object, the angle of view, the entrance 
pupil location, the relative aperture (f/number), the image 
distance, the shape of the focal surface (plane or curved), 
and (for some applications) the exit pupil location. The lens 
illustrated in Figure 1 has each of the prescription parameters 
labeled with its storage address number. Listed below are 
descriptions of the storage nomenclature and the conventions 
used as one starts at the first surface and proceeds through 
the lens. The lens prescription is stored in the memory by 
addressing the parameters as described below. The radii of 
curvature are sent to L132, L42, etc. Surfaces which appear 
convex aS seen from the left are positive; concave ones are 
negative; and the symbol for a plane is +0. The axial distances 
or thicknesses from the previous surfaces are sent to 181, L4l, 
etc. Distances proceeding to the right are positive and to the 


left, negative. The first distance (131) is from the entrance 
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pupil. With the pupil to the left the distance is positive, 
to the right it is negative. The aperture radii of the sur- 
faces are sent to L30, L40, etc. If the radii are not speci- 
fied the symbol +O is used, and the code will not check for 
ray vignettes from this cause. Central obstructions are spec- 
ified by negative radii. 

The optical surfaces are generated by conic sections. 

The eccentricities which describe the conic sections are sent 
to L33, 143, etc. (O for a sphere, 1 for a parabola, >1 for a 
hyperbola, between 1 and O or negative for axially thick and 
axially thin ellipses, respectively). 

The refractive indexes for the first wave length (color) 
are sent to L24 (object space index), 134 (first glass index), 
144, etc. For the second wave length the indexes go to L25, 
135, etc. Indexes for a total of six wave lengths may be used 
(to L29, L39, etc.). The refractive indexes may be used in 
any order and number desired by storing in L-89 the number of 
colors to be used, the number of the first color to be used in 
L-88, the second color in L-87, etc. (example: L-89 X3 X2 Xl 
x3). 

For a displaced and/or tilted spherical or conoid surface 
the aperture radius, the axial distance from the previous 
surface, the radius of curvature, and the eccentricity are 


stored at the usual addresses. However, some variation of 
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-88 
-83 


the normal procedure is necessary. Twenty storage addresses 
are required. In place of the first index of refraction, +0 is 
stored, a signal that more surface data are to follow. Plus 
one is stored in the next address, a signal that displacement 
and/or tilt parameters will be given. The next four addresses 
are cleared. The next 10 addresses have the x translation, 
the y translation, the sine of the angle of tilt in the x 
direction, the sine of the angle of tilt in the y direction, 
and the six indexes of refraction. (Example: 130 +0 +1.2 
44.65 +1 +0 +1 C4 +.015 -.01 +.03490 +.01745 +1.51462 +1.51700 
4+1.51901 +1.52264 +1.52712 +1.53047) 

A cylindrical surface (whose axis is parallel to the x 
axis and normal to the optical axis) requires a procedure 
which is similar to that specified for displaced and tilted 
spherical surfaces. The aperture radius, the axial distance 
from the previous surface, the radius of curvature, and the 
eccentricity are stored at the usual addresses. Twenty storage 
addresses are required. In place of the first index of refrac- 
tion, +0 is stored, a signal that more surface data are to 
follow. Plus two is stored at the next address, a signal that 
this surface is cylindrical. The next eight addresses are 
eleared, and the six indexes of refraction follow. (Example: 
130 +0 +1.2 +4.65 +1 +0 +2 C8 +1.51462 +1.51700 +1.51901 
4+1.52264 +1.52712 +1.53047) 


20 


SQqENTIFIC LIBRARY 


aug 27 1963 


U. S. PATENT OFFICE i 


-- 


ae 


- 


When a mirror surface is to be specified, a sign change 
is made in the refractive index of the medium which follows the 
mirror surface. The changed sign will be kept by succeeding 
refractive media until another mirror surface occurs, at which — 
time the sign of the refractive index is again changed. To 
specify surface radii and distances in mirror systems the same 


conventions are used as for refracting systems. 
RAY-TRACING PARAMETERS 


The distance to the object and the maximum angle of view 
are determined by the design specifications, or, in the absence 
of specific claims, the object distance is assumed to be in- 
finity (up to +1038), and the field angle is discovered by 
tracing series of ray bundles from object points at increasing 


angles. The object distance from the first lens surface is 


stored at L20. The angle of view is determined by the object 20 
height. The initial object height is stored at L21, the 21 
object-height increment is stored at L22, and the maximum 22 
allowable object height at L23. 23 
If the image distance is specified, it is stored at L14 14 
and the plane setting mode X1 stored at L-56. If the image -56 
field is curved, the additional sagitta for the first, second, 
etc., image points are stored at LO to L6. If the image dis- 10 
21 
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Si 


-49 
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tance is unknown, zero is stored at L14 and the plane setting 
mode XO is used. The ray-tracing program will then set the 
planes of interest at the back focus calculated from the merid- 
ional ray trace plus Li4 plus the corresponding sagitta in LO 
to 16. 

The distance of the entrance pupil from the first surface 
is stored at L131, negative if to the right and positive if to 
the left. If the entrance pupil distance is not specified it 
may be set at an estimated distance. 

The desired f#/number of the image beams is stored in the 
weight section at L-49. The symbol +0 is stored at L-49 if 
one wishes to fix the entrance-pupil radius when designing. 

If the f/number is known, the entrance pupil radius is hand 
computed and stored at L15. 

The numbers which designate the shape of the focal surface 
are stored in LO to Il6 and are zero if the focal surface is 
a plane. For a curved focal surface the numbers are positive 
or negative for convex or concave shapes, respectively. 

The image spot is evaluated from the size of the ray- 
bundle intersection on a plane normal to the optical axis. The 
predicted resolution is the mean deviation of the rays from 
their mean center. Any number of planes may be selected by 
Storing the number desired in L-57 (example: X3). The spacing 


between planes will be determined by the number stored in 
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Lll (example: +.01). The location of the first plane to be 


printed can be shifted by storing a number in L10 (example: -1 


sets the first plane one to the left of the plane of interest). 


The (paraxial) focal length is computed from a parallel 
ray traced automatically at an aperture height determined by 
the number stored in L1l2 (example: +.00001) multiplied by 
the maximum aperture stored in L115. 

The back focal length is computed from an axial-object 
ray traced automatically at an aperture height determined by 
the number stored in L13 (example: +.707) multiplied by the 
maximum aperture stored in L115. 

Feathering or negative thickness between adjacent lens 
surfaces is calculated with a paraxial ray and a marginal ray 
for the axial object position and with the upper and lower 
marginal rays for the extreme off-axis object position when 
the entrance pupil scale factor +1 is stored in L17 and a 
minimum feather thickness of +O is stored in L18. To avoid 
machine stops, wnich may occur if feathering develops, the 
scale factor and thickness values may be changed to smaller 
and negative amounts, respectively, at the convenience of 
the designer. 

The pattern of rays to be traced is determined by the 


numbers stored in L-170 through L-90. The pattern of rays is 


placed in a unit-radius entrance pupil. The ray-trace program 
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17 
18 


-170 


-169 
-168 


-167 
-~166 


Scales the unit-pupil dimensions to the actual pupil dimensions 
stored in L15 to obtain the coordinates of each ray to be 
traced. The rays are usually specified to sample the entrance 
pupil with a grid pattern in which each ray represents the 
surrounding area, but as many as 25 rays may be positioned 
individually in any pattern desired. Since rays traced through 
a centered optical system are symmetrical, only half the rays 
of symmetrical ray patterns need to be traced to obtain the 
ray intersections of the other half by a sign change of the 
x coordinate. Mode Xl stored in L-170 directs the program to 
trace the specified rays only and to generate the symmetrical 
ray intersections by a sign change in the xs (x = 0). Mode 
XO stored in L-170 directs the program to trace the specified 
rays with coordinates (x,y) and also rays with coordinates 
(-x,y) when uncentered optical systems are being considered 
(x is calculated). Mode X2 stored in L-170 directs the program 
to trace the specified rays only through centered or uncentered 
systems (x is calculated). The total number of rays to be 
traced in the ray pattern (100 maximum) is stored in L-169. 
The vertical ray spacing (Ay) is stored in L-168. 

The ray coordinates are now generated in columns by triads 
of numbers consisting of the x coordinate of the column of 
rays (stored in L-167, L-164, . . . L-92), the y coordinate 


of the first ray of the column of rays (stored in L-166, L-163, 
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L-91), and the number of rays in the column (stored in 
ielobmeb-16e, . . . L-9O). 

Figure 2 illustrates, and Table 1 gives the instructions 
for generating, the ray patterns which have been used for lens 
designing and for testing the performance of lenses. The two- 
ray pattern is especially useful during the preliminary design 
stage because it is economical of machine time while still 
evaluating adequately the major lens defects for most design 
problems. The six-ray pattern is usually adequate for the 
remainder of the designing, but a final run with the 12-ray 
pattern may occasionally greatly improve lens performance. 
For assessing lens performance, the many-ray patterns are 
especially useful when there is a widely variable vignetting 
of the beams which form the several images. 

The ray pattern on any specified plane may be graphed 
from the ray-intersection coordinates obtained from prints or 
punched cards which are supplied at the scale required by the 


point-plot factor stored in L16 (example: +100). 
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TABLE 1 
EXAMPLES OF INSTRUCTIONS FOR GENERATING SELECTED RAY PATTERNS 
Rays Card Image 
2 itatGoypal 42 +1 +.5 -.5 X2 
Two xl XS +1.2 +.35 -.6 X2 +.7 +0 Xl 
lO) xl Xee +1.4 +53 -.7 X2 +.% =3) xD 
+.7 +.3 Xl 
Pe iOuxele xS +.7 aes Sut 2G +.7 -.35 X2 
Gueel=awiO: Xi X6 +25 +.25 -.75 X4 +.75 -.25 x2 
Smet Oexde XG -peZt +.2 -.6 X4 +.6 -.6 X4 
lee EI70 X12 X1@ +.35 +.175 -.875 X6 +.525 -.525 x4 
+.875 -.175 xX2 
OMe Oo xa X16 +.3 +.15 -.75 X6 +.45 -.75 X6 
eG Bien tis) oh 
2Cm Or Xe X26° +225 +.125 -.875 x8 +.375 -.875 x8 
+.625 -.625 X6 +.875 -.375 x4 
VOMEE= 720k xXeLQ) 2 +.1 -.9 X10 +.3 -.9 X10 
+.5 -.7 x8 Stone soil 2K +.9 -.3 x4 
Cl OMX SGA 5 +.075 -.975 X14 +.225 -.975 x14 


L-161 


100 L-170 
L-161 
L-152 


+.375 -.825 X12 +.525 -.825 x12 
+.625 -.525 x8 +.975 -.225 x4 


Xl X100 +.125 +.0625-.9375X16 
+.3125-.9375X16 +.4375-.8125x14 
+.6875-.6875X12 +.8125-.4375x8 


+.675 -.675 X10 


+.1875-.9375X16 
+.5625-.6875x12 
+.9375-.3125Xx6 
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LENS-DESIGNING PARAMETERS 


The performance of a particular prescription is influenced 
by the position of the entrance pupil, the position of the exit 
pupil, the position of the focal surface, and the rays vignetted 
at the field edge. The optimum positions of the first three 
can be determined by the designing part of the code. Since 
ray vignettes from any cause will stop the designing sequence, 
the effect of vignettes can only be tested with the ray-tracing 
part of the code. 

To use the designing program it is necessary to provide 
(1) weights for the performance errors, (2) substitutions 
(if desired) for the lens parameters, (3) increments to the 


desired lens parameters, and (4) designing instructions. 


WEIGHTS 


The operating relative aperture (f/number) of the lens 


is stored in L-49. The entrance pupil radius will then be 


TUT wom 


recalculated before each design cycle. If the entrance pupil 


radius is to be fixed, zero is stored in L-49. The required 
principal focal length is stored in L-48. The weight on the 
focal length (if any) is stored in L-47. The exit-pupil posi- 


tion is calculated from a ray traced from an object point at 
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a height stored in L-45. The desired distance of the exit 
pupil from the last optical surface is stored in L-44. The 
weight on the desired exit-pupil distance is stored in L-43. 

Optimum performance generally requires minimum image-spot 
sizes on the plane of interest. These image spots -are made 
up of a maximum of six colors each from a maximum of seven 
object points. The weights for these spots are stored in L-39 
to L-33. Weights on the lateral chromatic aberration for each 
of these seven images are stored in L-29 to L-23. Weights on 
the desired image height for the average height of each of the 
seven images are stored in L-19 to L-13. The seven desired 
image heights are stored in L-9 to L-3. The above weights 
are the total number which control the design program. The 
number stored in L-46 determines another control factor, the 
maximum parameter changes allowed by the program. To obtain 
the individual limits, this number is multiplied by each 


parameter alteration. 
SUBSTITUTIONS 
The code permits the automatic designing of a lens in as 
many as five conditions of use (for example a zoom lens) by 


means of a substitution routine. The addresses (L-200 to 


L1499) of the parameters to be substituted are stored in L1500 
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1500 


1600 


1010 


to L1599 (example: L1500 X14 X-49 X51). The corresponding 
parameters to be substituted first are stored in L1600 to 
L1699 (example: L1600 +350 +6 +100). The second set of 
parameters is stored in L1700 to L1799 (example: 11700 +400 
+7 +150). The third set is stored in L1800 to L1899 and the 
fourth set in L1900 to L1999. The substitutions will be made 
and designing accomplished automatically if the number stored 
in L-59 is the number of sets of substitution parameters and 
the number stored in L-58 is the number of substitution param- 


eters in each set. (example: L-59 X2 X3). 


INCREMENTS 


Any of the prescription parameters (LO to L1999) may be 
altered by the design program to give optimum lens performance. 
Storage has been provided for the design of as many as 49 
individual parameters. As many as six parameters can be given 
identical alterations simultaneously in any plus or minus 
combination. The increment instruction requires 10 storage 
cells, starting at L1010, L1020, ... LI490. The first num- 


ber, X1 to X6 or X-1 to X-6, is the number of addresses on 


mt i a 


which this increment is used. The sign of this number is a 


convention. If it is positive, the increment is added to the 


parameter; if negative, the increment is added to the recip- 
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rocal of the parameter (for radius of curvature). The second 
number, XO to X1999 or X-O to X-1999, is the positive address 
only of the first parameter which uses this increment. The 
sign of this and the next five numbers is a convention. If 

it is positive, the increment is added to the parameter or its 
reciprocal; if negative, the increment is subtracted from the 
parameter or its reciprocal. The third number, XO to X1999 
or X-O to X-1999, is the address of the second parameter which 
uses this increment. The fourth, fifth, sixth, and seventh 
numbers are the addresses of the third, fourth, fifth, and 
sixth parameters which use this increment. All parameters 
addressed here will be given the same amount of alteration, a 
capability useful for the design of symmetrical systems and 
also for holding selected lengths or curvatures equal or con- 
stant as desired. The number stored in the next cell is the 
program-computed single determinant for the increment used. 
The next number is the single determinant specified for this 
increment by the designer. The last number is the actual 
increment used by the design program. The minimum value for 
any increment (example: .000001) is stored at L19 and is 
automatically transferred to L1019, L1029, etc., if the incre- 
ment calculated from the required determinant is less than 
the assigned minimum. Five examples of parameter-inerement 


instruction card images are shown in Table 2. 


Syl 


1011 


1012 


TABLE 2 


EXAMPLES OF PARAMETER-INCREMENT INSTRUCTIONS 


TOMO) Xa) X82 LIOI8 +.12 +.0001 

imiolOmex—d XSe C6 +.1 +.0001 

L1020 Xe KOM aC St. Ol “Ee 1 

L1030 X-4 x42 x52 X-62 X-72 C3 +.1 +.0001 
Li0#0 X-6 X32 XHe@ x52 X-62 X-72 X-82 Cl +.1 +.00001 


The increment-instruction sequences assume the serial identi- 
fications 1, 2, 3, etc., which are used to call any one or 


any continuous series of them for designing. 


DESIGNING INSTRUCTIONS 


The design-parameter increment-instruction sequences 
described above can be used only if the required instructions 
are stored in L-79 to L-65. The minimum identification number 
to be used is stored in L-79 and the maximum in L-78. The 
number stored in L-77 is the number of times this set of in- 
erement cycles is to be used for designing. The total number 
(1 to 10 permissible) of parameter increment identifications 
to be used in each cycle is stored in L-76. The sequence of 
identification numbers to be used is next stored in L-75 to 
L-66. With each design cycle each number in the calling 


sequence is incremented by 1 (any becoming larger than the 
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maximum will be reset to the minimum) so that all parameters 
in a selected series will be used in serial order and repeated 
if enough design cycles are ordered. This arrangement is used 
because it is frequently not possible or desirable to use all 
parameters at one time, yet all must eventually be allowed 

to work the prescription. It has been found advantageous to 
arrange the increment calling sequence so that one obtains 

the greatest variety of interactions between parameters. With 
crude starting prescriptions it is often found that the single 
and multiple determinants tend to fluctuate violently, and for 
this reason designing is usually started with only one or two 
increments per cycle. As the design is stabilized the number 
of increments per cycle is increased. Most designing has been 


2 to 10°- 


done with the single-increment determinant in the 10 
range. With eight increments being used simultaneously, the 
determinant frequently ranges from 107? to ior =? When the 
determinant is in the Ones to 10739 range, there is often 
great fluctuation in the design progress. If a particular 
combination of increments does not improve the prescription, 
the use of that combination will be skipped for the same num- 


ber of cycles as the number stored in L-65. Six examples of 


design-cycle instruction card images are shown in Table 3. 
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TABLE 3 


EXAMPLES OF DESIGN-CYCLE INSTRUCTIONS 


L-79 X1 x12 X25 X1 X1 C9 

fevomet x42. X37 X2 X1 X5 c8 

Teqoma xe x18 X4 Xi Xe x4 x8 c6 

Reomedeexdon X10 X6 Xl X2 X3 X5 X7 X10 c4 

yom x29 X6 X1 X2 X3 X4 X8 X11 ch 
PevCmaueaaxeo XS X5 X6 X7 XO X12 X13 X14 xe ce XI 


PARAMETER INTERACTIONS 


Because of the interactions among parameters, the design- 
progress rate increases with the number of increments working 
simultaneously. The number of interactions between n param- 
eters working simultaneously is 3n(n-1). (For example, with 
this program one obtains 0, 1, 3, 6, 10, 15, 21, 28, 36, or 45 
parameter interactions depending on the number of increments 
worked simultaneously.) The design-progress rate is increased 
as the number and variety of the parameter interactions in- 
erease. The working (incremented) parameters are serially 
stored, and the calling order is at the command of the designer. 
For example, with a series of 15 parameters in which five are 
being worked simultaneously, one might select every third 
parameter, and the three calling sequences for designing would 


PPro omlOmess 2, 5,0, 11, 14; and 3, 6, 9, 12, ‘15 yamAun 
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parameters are used every three design cycles, but only two 

varieties of parameter interaction are present, i.e., those 

spaced 3 and 6 apart. At the other extreme one can select the 

M5vcalling sequences: 1, 2, 3) 4) \5:) sc enGMel > pile OREEeE 

In this case all parameters are worked every 11 design cycles, 

with four varieties of parameter interaction present. Or one 

might use the sequences 1, 2, 3, 5, 14; 2, 3, 4, 6, 15; 

Here all parameters are worked every nine design cycles, with 

five varieties of parameter interaction present. Or the re- 

verse sequence 1, 3, 12, 14, 15; . . . would be equally useful. 

However, the greatest variation of interaction is obtained 

with the 15 calling sequences 1, 2, 5, 10, 12; 2, 3, 6, 11, 

13; . . . , with all parameters worked every five design cycles, 

and all seven possible varieties of parameter interaction 

present. For clarity this information is summarized in Table 4. 
TABLE 4 


ANALYSIS OF PARAMETER INTERACTION IN SELECTED SEQUENCES 


Parameter Spacing 
Incidence Frequency 


Spaces Between 


Desi 
Sign Parameter Parameters 


Calling Sequence 
is LOseals 
1, 2, 35 4, 5 

9 27 35 Sy ae 


Fr I nO = Oe 
pet (Gay KU) | (U8) 7S) 
fee) Se SL 
TOS EDO RAD hess Oe 
(i), Gia ey Qe 
pep fey SH 
W Cy “a Cy ene 


al 
ty 3, tes Wh, WG 
1 


Dp LO), je 


0 


Table 5 gives an extensive series of parameter calling 
sequences, which have been selected in each case to provide 
all possible parameter interactions. The number of interac- 
tions of each type is distributed as evenly as possible, and 
each parameter is repeated as often as possible by avoiding 


sequences with large parameter spacings if there is a choice. 


SAMPLE CALCULATION 


The lens to be designed in the following sample calculation 
is similar to the lens illustrated in Figure 1. It consists 
of two cemented doublets in series which produce an F/2 Lister- 


fi 


type lens’ set to work at a magnification of about one third. 
(To increase the field of view and field curvature for this 
illustration, the space between the two components was reduced. ) 
The objects are at an axial distance of 12 inches from the 
first lens surface and 0, 1, 2, and 3 inches from the optical 
axis. The images are on a curved field which is held at an 
axial back focal distance of 3 inches by the use of plane 
setting mode Xl. The design was started with two cemented 
sheets of 651584 and 720362 glass spaced 1 inch apart. The 
entrance pupil distance was fixed at -1l inch and the aperture 
radius at 1 inch. Figure 3 shows the complete series of card 


images used to design this lens. The data for the starting 


prescription, the weights, the ray patterns, the parameter 
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L- 1948 NSSAPPLE PROBLEY FOX ILLUSTRATION s 

L-19? HS O1--30--63 5 

L-89 XS KI X2 XB) U-57 Xe KT L-86 #20 L-39 #100 #100 *100 +100 

LIC -2 #23 e.0001 +2707 #3 ©) ©1000 #1 -1 #.N001 

L2G 22 oO eh eS.1 2 eb eh LSU eO -t 60 #0 #16651 e1eouTTT 1.0508 
LOO 00 1.25 oC oO e1.7T? 1. TIS 1.73422 150 40 #25 #0 00 #l ele 
LOO #0 ef #0 #0 1.651 *1eoerTh #1265885 

LTO 60 1.25 60 00 e172 e1.7IRSH C1L73B272? LBO 40 #25 40 40 el el & 
LiClO x-1 X52 CO ’ OY £1020 X-1 KS2 CO #21 *.ul K-1 X67 CO eet 
L10ao x-) x52 C4 Ol xX-1 Xb? CH #.) #202 K-1 X72 CO #.d #201 
LICTO KT XE CO 2.05 #21 KI X2 CH 0.05 O21 XE KS CH 6.05 Oe) 

LOO KI x31 CO e205 oT 

L-170 KTX2 #1 0.5 =25 KZ U-79 XI Xe xXy XI KI CY *2 
a-7¥ KAD XO KIS) ato KI CY ? 

c-79 XY xv xIG@ xr xT otY Lyd 

U8? ee L-79 XI XY x19 X2 XI XN CU *2 

t-79 AL XY X28 xh XI x2 X5 X7 «CH +2 

L-170 XI KO #65 #225 -275 KW 6.75 =.25 K2 L-TYK VAY LUKSKIX2KSXSCO 
L-79 KU KIO RED XO XT -X2 KS XB rH XK Ch 82 

L-U7O X1 X12 +.5555 +.1657 -.9555 KO #25 ~25 X& %.6333 -. 1067 X2 
L-t¥ KP KIO XS XH KI X2 XS XB XO XY Ch 2 
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SAMPLE PRORLEM 


7.0699999-01 
1.0000 
3.1000000+00 
0.0000U0U+00 
9.0000000+09 
0.0000U00+0U 
0.0000000+00 
0.0900000+0uU 
0.0000000+00 
SENSE SRIIC 


2.0000000"u1 
§.0U00000¢02 
0.00000000u9 
0.0000000+L0 


FOR ILLUSTRATION 


3.000000u+0U 
V00+00 
1.00G000u+00 
1.651000U+00 
1. /2u0000+00 
1.90000000+00 
1.6510000+00 
1.7200000+00 
1.000000U+00 
HH MODE 7ERO. 


0.00U0000+00 
0.G0000U0+00 
0.0U40000+ v0 
0.00000U0+90 


~1.9000000+00 


PRESCRIPTION 
19 7220600000900 ¥.u000000-01 9. 99999T9-05 
20 1.2000009+C) 9.0000000+C0 1.00000UU+00 
30 0.0000000+CO =1.0090000+C0 4.909000U0+00 
eo 0.v0U9G00"00 1.2590000+C0 0.0U00U00900 
50 0.0000000¢50 $.u0000U0-C1 0.000000 
‘ oo 0.v000000+00 1.0000000¢C0 0.060090 
_ 1 0.0000000+00 1.2500000¢00 0.0000000°00 
su 0.0000000+00 5.0000000-C1 0.0000000¢CG 
PLANE SETTING MODE Ove. LATIECE MOOE OuE. 
wElonTs 
woh O.0G000C09L0 0.9000C00+C0 0.0000000+00 
é -s¥ 1.00000L0¢02 1.0000000+C2 1.00U0000+02 
“29 0.0000000eC9 0.0000000+00 O.UL00U0090u 
19 0.0000NC0+CO 0.0000C00+CU 0.00000U0+0U 
-¥ 0.u0000u0+00 ¥.0000000+00 0.0040000+09 
i LATTICE N= 2 DELTA Y = 1.0000000+00 
$.0009000-01 -5.0900009-061 2 
INCREMENT IDENTIFICATION... 1 ue % 1 
COLOR IDENTIFICATION... nf 9 (ay >} 


0.0000000+u00 


¥.0000000400 


*2 


1.0000000+G0 


1.0000C00+00 
1.64/7100+00 
127145400400 
1.0000000+00 
1.6477100+00 
1.714 5400"00 
1.0000UU0+00 


NUMBER OF PLANES 


0.GOv0000+CU 
0.v000000+u00 
0.C000000+00 
0.C000000+00 
0.cv00000+00 


O1--30--63 


1.0000000+¢03 


9APPPPVTS“O5 


1.0000000+C0 
1.6588500+00 
1.7542200+00 
1.0000000+00 
1.6588500+00 
1.7342200+C0 
1.0000000+Cc0 
4 


0.c00C000+90 
0.G000U00¢0U 
0.0000000+90 
0.0000000+00 
0.00u00000+00 


10. toc. NO. ADDRESSES OF DESIGN PARAMETERS ACTUAL DET. RECUIKED LET. INCREMENT 
1 1o1o =7 32 G.900000uU+00 1.00000C0-0)1 a VIVIIGI-O$ 
2 1020 | 52 0.0u00G0U+00 1.00000C0-01 FaPPYIVII-U% 
3 1030 <1) 02 0.v0U0U0U0+00 1.000U0C0-01 FeVIVIIIV=-08 
“ Toso -1 62 0.000000u+00 1.0000000~u1 2. 99499999=05 
f 5 1059 -1 42 0.000000u+00 1.00000C0-01 YeIPIIVII=O4 
i 6 1060 =i T2 0.n0U000UFCU 1.0000000-01 ¥eVIVIPVI-OS 
f oro 1 1 ©.Y00U0U0900 = 5. VOLO000-v2 1.0C00000-6) 
4 1080 1 2 U.J00000U+0U 5.0000000-07 1.00000UU-u01 
9 1O?u 1 3 0.000000L+00 5.0000000-02 1.0000000-01 
$ 1) 1100 1 31 0.000000U00 5.0000000-02 1.0000000-0) 
Ceres AAS a6 58! (9) 20 INCK. FACTOR  DETEKMIWANT OLD LENGTH NER LENGIN IMPROVEMENT 
mEXIT 1 5.4053d47-01 A.N56B02600) 3.403 (419002 1.5550818+02 5-4269074-01 
MERIT 2 1.0090000+00 6.5409 19S¢01 1.555d818+02 TeVNN7392¢0R 9%. 2854603-01 
Mi mrall 3 1.0000000+00 2-1ATI9SS*O V.EDT7392e0R 929965106900 1.0062233-91 
meat 6 1.000000U+06 T.AVOVIS 100 ¥.¥965190+00 9-9636626+00 1.2852504-03 
MERIT ' 1.00000U0+00 8.7207531-02 929836626400 9.9557289+U0 2.09794 26-03 
BERIT 2 1.0000000+00 1.1838008-01 ¥.¥5577284+00 9545368400 1. 1975939-0% 
meet 2 1.0000000+00 1,003217Ty¥-01 ¥.¥5S4536U400 9. 9462004460 6. 374K 979-04 
MERIT « 1.0000000+00 1.0070031-01 929462004400 9.9337771+00 1.2490475-03 
MERIT ! 1.0000UU0+00 9.7SvTe2b-02 ¥-93S37TTTIVOO 9.90586012000 2.8101934-03 
PRESCRIPTION SAMPLE PROBLEM FORK ILLUSTRATIUN 01--30--63 
a0 ~2.0000000+00 1.00000U0-01 %e99PTITS-“O5 7.0699999-01 3.0000000+00 1.0000000+00 1.0000000+03 
1.00U0000+00 ~1.0000000+00 9 AP99V9TS-05 
ae 1.2000000¢01 yoouuy0+00 1.0000000+00 5. 1000000+00 1,000000u0+00 1.0000u00+00 1.G000000+C0 
30 G.0000000900 = -1.4u000U00+00 5. OK3847TH400 0.0000U00+00 1.05100UU+00 1.6477100¢00 1.6586500+00 
ao 0.0900000+00 1.¢500000+00 G.0000U00e00 0.000000U+0U 1.720000U0+00 1.7143400+00 1.7342290+C0 
4 ou O.vuou00I¢L0 5.vu00000-01 922706752400 0.0900000+00 1.U00000uU+0U 1.0000000+00 1.0000000+CU 
“od 0.0000000+0u 1.009000u+CO ~1.O48THES*U2 0.U000400+00 1.45100UuU+00 1.6477100+00 1.6588500+CO 
Led 0.0000000+0U 1.2590000+Co 0.0000uLu9L0 9.9000000+90 3.7200000+00 1.7143400+00 1.7342200+00 
ov 0.0000690+00 5.G000000-01 ~2.82VK 266002 0.u900U00+00 1.uouGuUu90u 1.0000000¢U0 1.9000000+00 


Figure 3. 


Printed output 


showing Lister lens design start. 
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alterations, the design instructions, and the machine opera- 
tions are given. This information is followed by the complete 
prescription data print, the accomplishments of each of the 
next nine design cycles, and the resulting lens prescription. 
Designing was started with three colors (C, D, F), a two- 
ray pattern, and a series of nine single increments to the four 
glass-air surfaces. As a result the prescription's figure of 
merit went from 340 to 9.9. A summary of the complete 9-minute 


series of runs used to design this lens is shown in Table 6. 


TABLE 6 
SUMMARY OF LISTER LENS DESIGNING SEQUENCE AND PROGRESS 
Operation Ray Increments Used Cycles Figure 
Performed Pattern “Total Per Cycle Run of Merit 
Start 2 2 2 - 340. 
Designing 2 4 1 9 9.9 
Designing 2 6 1 13 9.7 
Designing 2 fe) 1 19 5.9 
f/2.0 2 - - - 15.4 
Designing 2 9 2 19 6.6 
Designing 2 9 4 28 16 
Ray Change 6 - - = eG 
Designing 6 9g 4 10 2 
Designing 6 10 6 11 0.96 
Ray Change ie - - = 0.96 
Designing 12 10 6 5 0.94 
Se ido J 
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These runs produced an f/2.0 lens of 3.461-inch (D-light) 
focal length, which should give 400-lines-per-inch resolution 
on a 3.0525-inch-radius concave focal surface. There is a 
maximum of 0.002 inch lateral chromatic aberration at the edge 
of the 13-degree half field of view. Figure 4 shows the print 
of the final prescription and the 12-ray analysis of lens 
performance. The spot diagrams of the four images are shown 


in Figure 5. 


A COMPLEX DESIGN 


An example of a complex design which can be made with 
this program is a special-purpose zoom lens. This lens is an 
objective for a high-speed camera used to photograph explosions 
through the 4-inch-diameter viewing port of a concrete bunker. 
Twenty-millimeter images of 150- to 450-mm explosive objects 
are required, with the objects 9000 mm from the front lens 
element. To obtain this range of magnification, a zoom ratio 
of approximately 3.6:1 is needed. Between the viewing port 
and the objective's image, there is a distance of more than 
2000 mm, which can be used for any desired arrangement of lens 
elements. 

in the search for a solution to this difficult imaging 


problem, the lens designs of US Patents 2,778,272, 3,051,052), 
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PRESCRIPTION FROM DUMP NUMBER ! 


SAMPLE PACBLEM FOR ILLUSTRATION 


O1--30--63 


9.0000000e09 G. 000006000 


0.0000000+00 = ~1.7951207-02 — -6.6173003-02  -1.42430/2-01 6.006000u+00 o 
10 -2.0000000+00 1.u000000-01 229999975-U5 7.0699999-01 3.0GU000N+00 160180819005 2. Co900G0854 
1.9UU9000+00 -1.0900000+00 = ¥. 99 99975-95 
20 1.2000U00+01 0.0000000+00 1.0000000+00 3. 1000000+00 1.G000G00+00  1.0000ug0+00 1. 0u0ga90009 
30 0.000U000+00 -1.0133081+00 6-9734230+00 0.0000000+06 1.6510000000 16477100605 «4. 6SBBSOD4 00 
nO 0.¥000000+00 1.250000U+00 = - 1.5508849+00 0.G000000+09 1.7200000+00 =. 185890000 ©=— 1 7 H¥2ZOO8CS 
50 0.0000000+00 5.V0U0000-C1 -6.9003449+00 0.G900000+00 1.090000u900 =: 1. 0600500008 
00 0.0000000+00 1.0000000+00 8.2007210+00 0.0000000+00 1.6510000400 = J.44877100800 1658500009 
70 0.0U00U0U+00 1.2500000¢00 =~ 14668796400 0.v000000+00 1,720000U*0G =. 7143400000 «= PNZZOO4CD 
30 0.0000U00+00 5.0000U00-01  -3.6679682+0U 0.0000000+0u 1.0000000+CO = -:*1.0009u00%00 ©=—- 1. Ca99000+G2 
PLANE SETTING MODE ONE. LATTICE MODE ONE. SENSE SWITCH MOCE ZERO. NUPBER OF PLANES 
CoLoR FOCAL LENGTH FOCAL POINT H EXIT Puri E/NUMBER BACK FOCUS RAYS 
1 34609916400 1.9274 388+00 0.000000uU+00 2.0006060+00 2.99343 TEVOD 12 
LOC UF PLANE AVERAGE X AVERAGE ¥ RYS K eas ¥ SPST Sile 
2.4000000+00 0.0000000+00 -2.3685035-06 2. 2382708-02 2.2330089-02 36165287192 
2.7000000+00 Q.U000000+00 ~1.10%4321-06 1.085080s-02 1.GbS 9812-02 1.5388 C37-02 
3.0000000+00 0.0000000+cu 1.5964473-07 1.3199177-03 1.3177530-¢3 1..3665503-93 
3. 1000000+00 0.0000000+00 1.4 237206-06 1.237005 1-02 1.2368903-02 1.7893082-02 
2 324651695400 1.9255137+0C 0.0000000+00 2.0094%4585+00 229896925000 12 
2.80v0000+00 0.9000000+00 -2.2860038-06 2.1875057-02 2.1873056-02 3.0938595-02 
2. 2000000400 0.0000090+00 -1.0217773-06 1.0356273-02 1.9355341-02 1.8545323-32 
3.00000u00+00 0.0000000+u0 2.4245432-07 1.7874692-03 1.7872633-03 2.5277 175-93 
3.1000000+00 0.9000000+00 1.5066859-06 1.2877522-02 1,2596305-u2 1.8233990-02 
4 3.459 7524 +00 1.9371255+06 0.0000000+00 2.0024781+00 B.GNloa72°00 12 
2.HU0UV009+00 0.0000000+co -2.1098770-06 2.4650463-02 2.4653501-02 3.827%057-92 
2.9000009+00 0.0000000+¢c0 1.449 259-06 1,3125655-02 1.$124356-02 1.2561561-02 
3.0000U00+00 0.u000000+00 -1.8399036-07 1.6848 162-03 146885927-C3 2. $225 519-95 
3. 1000000+00 0.0000000+v0 1,0718672-06 9.99427 12-03 9. 9908499-03 ¥.0129211-02 
1 1.000u000+00 3. 6064031400 12 
2.7820485+00 v.0000000+C0 3.0654871-01 242792981-02 2.2360278-02 5.1929 681-02 
2.8820643+00 0.00000GU+U0 3.1158678-01 1. 1232371-02 1.0770105-02 1.9979425-92 
249820487400 0.0000000+uu0 3.1622576-01 1. 1073005-03 1.1907757-03 1,.62603975-O5 
3.0820487+00 0.0v00G00+00 3.2106 354-01 1.2029262-02 1,24%8COS-C2 1, 7310576-02 
2 1.0000000+00 =3. 6135179400 12 
2.78204 48+00 0.0000000+00 3.0640222-01 2.2275349-02 2. 1820980-02 3. 11E2849=02 
2.882U48n+0U 0.v000000+c0 3.1143812-01 1.07276 79-U2 1,026 14s4-c2 1, 484% 795-02 
224820447400 0.0000C00+c0 3.16474U2-01 1,5204409-03 1.7 1601TO-C3 2.2926 969-05 
3.0820487+00 0.0u000u0u0+C0 3.2130992-01 1.2566453-02 1.3003705-02 1, 8085475-02 
3 1.0000000+00 -3.5971508+00 2 
2. 18204858400 0.00v0G00+CU 3.0602654-01 2.5098595-02 2.4 719955-02 3.5225201-02 
2.d5204H48+00 0.0000000+00 3.1087417-01 1. 35533551-02 1.3154e68-C2 1.8373306-2 
2.¥8204687+00 0.0000000+00 3.1572180-01 2.U354%25-03 1.6469301-C4 2.5120045-25 
3.UH204UK7+CU 0.0000U00+CU 3.2056943-01 9.0279 156-03 9.9990 360-035 1, $859351-02 
; 2.0004000+00 =3.6225702+00 12 
2.7338270+00 0.v000uu0+Cc0 6.0558410-01 2.349% 662-02 2.1691736-02 3.211303e-02 
2..8338270+00 0. 0000C00+Co 1518953-01 1. 18545 18-02 1.02235%2-02 LSE 
2.9338270+00 0.0000C00+00 6.2479696-01 9. 377590?-us 1.536070I-03 1. 79VO858-03 
3.0335269+00 0.0000000+00 6.3440039-01 1, 1533753-02 1.3132303-G2 1. 187312802 
o 2.0000000+00 -3.6297793+00 12 
2.7338270+00 0.0000000+00 6.0609328-01 2.27876 10-02 2.1295873-02 3. 1279318-02 
2.4338270+0U 9.0000U00+¢60 6.1569399-01 1.131757 1-u2 9. 5V829 30-05 N.4a37Car-92 
?.9334270+00 0.0000000+00 6.2529409-01 1. 1650 930-03 2.1590160-05 2.4533226-03 
$.0333269+00 0.0000C00+cu 6.3489580-01 1,2099325-02 1.377 1157-02 Ve8331351-02 
= 2.0000000+00 -3.603593U000 12 
- ; aa -09 
2.7336270+00 0.0000000+60 6.U453108-01 2.58679 17-G2 2.8500787-02 3.825488 
2.0338270+00 0.0000000+CO A.NLISS31-01 1.42419 19-G2 1. 2853325-02 Lda 
2-9338270+00 0.0000C00+00 6.2377 770-0) 2.6287653-03 1.4 107292-03 Zoe laoee 
3.0336269+00 0.0000400+C0 6.3340407=01 2. 00Svaly-03 1.0631516-02 1.3853 456-02 
, 3.0Uv000G00Y = 36511435020 12 
246575693400 0. G0uUCU0+C0 #9051340 244571533-02 2.1629125-02 3.2605222-02 
247575093400 0.000G000+CO %.OKPKIGI-O1 1.2810 787-02 9. 1058506-05 Lease 
2.45750924C0 0.0U09UG0+00 ¥.1dv8 144201 1.31¥1555-03 2.7113205-C3 peo laats oy 
24575692400 0. 0000C00+00 9.3321609=01 HeuT¥71ba-02 1. 8OP4 380-02 1. 2TSS929~ 
e %.uuououg+an ~3.65380822+00 12 
2.0575693+00 0.U0v0000+c0 d.9128233-01 2..3975023-02 2.0627262-02 Rosettareras 
2-75756934C0 v.0000C00+C0 9.055095 1=U1 1, 2221248-02 G.d92 1810-03 DeSTESaT I-02 
2.d5756¥2+CU 0.0000000+00 V.LYT3009-01 1.0780 785-03 5. $028040-03 34 TAIOI0-05 
2-4975692+6u 0.0000C00+00 923390347-01 1.14) 1086-02 1, WROOUSe-02 1.3 laQo?3-02 
=~ 3.00000NG+u0 ~3.0318620+00 2 
226575093400 0.0000C00+00 .5591515-01 2..7062503-02 2.8460 13T~02 3.0881 133-02 
2-7575693+00 0.vv00u00+cu 9.0317 748-0) 1,5313308-02 1.2789176-U2 1, 89aVS20-U2 
228515692 +CU 0.0000C00+c0 Yah (U3973-O1 3.5745015-03 2.3193620-03 8.2778865-03 
¥23170202-01 Bo 2249583-03 DVM ITI3-92 1, 3850029-02 


2-9575092+C0 


Figure 4, 
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Printed output showing Lister lens performance. 
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Figure 5. Spot diagrams of four images made by Lister lens. 
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and 3,057,259 were tested and all were found to be inadequate 
for our purposes because of insufficient correction of off-axis 
aberrations and/or excessive shift of the image over the zoom 
range. It was then decided to let the lens-designing program 
find the optimum design for lens systems with three or four 
or more elements so that we could select the simplest design 
which would give the required resolution. Using the program's 
Substitution capability, we let the machine design each lens 
system simultaneously for five equally spaced positions of 
the moving elements. Several monochromatic designs were made 
with zoom arrangements which previous workers had found to be 
useful. The so-called four-lens system, which has five ele- 


ments, was found to be the minimum which would give the desired 


correction of all aberrations. The entire available length 


was used because a long lens nearly always gives the best 


performance. The design was then restarted with five cemented 


sheets of 517645 and 617366 glass (C, e, and g indexes) in a 
zoom system consisting of two coupled sliding lenses which 


move between three fixed lenses. ‘These paired components 


quickly assumed shapes to give negative, positive, negative, 


Positive, and positive powers. Before the final design was 


achieved, several readjustments of the movements and spacings 


of the elements were necessary. 


The prescription data print is shown in Figure 6. Table 7 
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PRESCRIPTION FROM 


DUMP NUMBER 53 


BRIXNER 700M FOR &/1 POSITION ONE. 


01--31--63 


Yo ~1,0000000+00 2. 999V979-03 1.90000000-05 7.0699999-01 3.76535234+02 5. 4UUS9AIeGl 1..0000000+00 
1.0000000+00 —2.0000000+01 $.9999999-0T 
Pad) 9.0000U00+C3 9.v000000+00 3.7500000¢G1 7.6000000+01 1.9002060+00 1.0002060+CO 1.0002060+00 
1.0002070+00 1.0002080+00 1.0002100+00 
30 0.0000000+00 -8.C000900¢02 5.0864065402 9.0000000+00 12515210000 1.5177200+00 1.5196600+CO0 
1.5233500+00 1.5276900+00 1.5512600+00 
&o 0.0000000+S0 1.0000000+01 2.95856ubr02 0.0000000+00 1.6117600+00 1.6167200+00 1.6206900+00 
1.6286000+00 1.6363900+00 1.6469200+00 
so 9.0000000+00 1.5000000+0) 3.50U59807+02 0.0000000+00 1.u00206u+00 1.0002060+C0 1.0002060+00 
1.0002070+00 1.0002080+00 1.000210u0+00 
60 0.0000000+60 1.U000000"C1 3.5075321+62 0.0000000+00 1.5152100+0U 1.5177200+00 125196600400 
1.5233500+00 1.5276900+00 1.5312600+00 
70 9.0000000+00 2.0000000+01 =6.6271720+02 0 .0000000+00 1.61317600+00 16167200400 1.6206900+00 
1.62860000+00 1.6383900+00 1.6469200+00 
30 9.v000000+CU 1.9000000+01 5.9126555+03 0.0000000+00 1.0002060+00 1.0002G60+U0 1.0002060+00 
1.0002070+00 1.0002080+00 1.0002100+00 
90 0.0000000+C9 5.u000000+02 625545400404 0.0000000+06 1.5152100+00 1.5177200+00 1.5196600+00 
1.5233500+00 1.5276900+00 1.5312600+00 
1090 0.0900000+00 1.0900000+01 1. 1846093402 0.0000000+09 1.6117600+00 126167200460 1.62069700+00 
1.6256000+60 1.6365900+00 1.6469200+00 
lio 0.v000000+00 1.5000000+01 2.1867210+02 0.0000000+00 1.0002060+00 1.0002060+00 1.0002060+00 
1.0002070+00 1.0002080+00 1.0002100+00 
120 0.0000000+00 1.0000000+0! 5.8100050+02 0.0000000+00 1.5152100+00 1.5177200+C0 1.5196600+00 
1.5233500+00 1.5274900+00 1.5312600+00 
130 0.00900000+00 2.v000000e01 =2.2635b0T9+02 0.0000000+00 1.6117600+00 1.6167200+00 1.6206900+C0 
1.6286000+00 1.6385900+00 1.6469200+00 
seo 0.0000000+09 1.40C0000+01 -6.78B35659*02 0.0000000+00 1.0002060+00 1.0002u60+U0 1..0002060+00 
1.0002070+U0 1.0002080+00 1.0002100+00 
159 0.0000009+00 1.6000000+03 &.b93TSI3+02 0.0000000+00 1.5152100+00 1.5177200+00 1.5196600+00 
1.5233500+00 1.5276900+00 1.5312600+00 
160 0.0000090+90 2.u000000+01 =2.2168650+02 0.0000090+00 1.6117600+00 1.6167200+00 1.6206900+00 
1.6286000+00 1.63283900+00 1.6469200+00 
170 0.0000000+99 3.yo00000+C! -6.4920066+02 0.0000000+00 1.0002060+00 1.0002U60+00 1.0002060+00 
1.0002070+00 1.0002080+00 1.000210C0+00 
PLANE SETTING MODE OnE. LATTICE MODE ONE. SENSE SWITCH MODE ZERO. NUMBER GF PLANES 3 
WEIGHTS 
~4y 1.20N0000¢C1 #.0000000+02 0.0000000+00 1.0000000+G2 0.00u0000+C0 0.co000000+00 0.0000000+00 
-39 1.0U90000+02 1.0000000+U2 1.0000009+02 0.0000000+v0 0.U0U0000+0u 0.0000000+00 0.0000000+00 
-27 0.us099000+00 5.9000000+C1 1.0000000+02 0.0000000+V09 0.00u0000+09 0.0000000+00 0.0000000+00 
19. 0.0000090+00 ©.0000U00+LO 1.0000000+02 0.0000000+Uu 0.0000000+00 0.0000U00+UU 0.U060000+00 
= 0.yoou0c0+00 0.0000000+00 1.0000000+01 0.0000000+00 0.0000000+00 0.0000000+00 0.0000000+00 
LATTICE N= 6 DELTA Y = 5.0000000-01 
2.9000000-01 -7.5000000-01 u 
7.500000u-01 -2-590C000-61 2 
INCREMENT IDENTIFICATICN. ~~ 1 18 12 10 8 yo afop RPS 2!) AR ice 4 
COLOR 1DENTIFICATION... 3 3 5 
10. Loc. NO. ADDRESSES OF DESIGN PARAMETERS ACTUAL DET. REQUIRED DET- INCREMENT 
1 1010 -1 32 1.0001991*02 1.0000000+02 1.1339532-06 
2 1020 -1 42 1.9000682+02 1.0000000+02 45..7589209-06 
3 1050 -1 52 1.0465263+02 1.00000U0+02 9.9999999-07 
‘ 1049 -1 62 1. 0000476+U2 1.0000000+02 1.0647986-06 
Ss 1090 mil 12 1.0000028+02 1.00000U0+02 5.4958878-06 
6 1060 = 62 1. 1893105402 1.0000009+02 9.9999999-07 
7 1076 = 92 1.0000299+02 1.0000000+02 2.0959069-06 
e 1030 -1 102 1.0000603+02 1.0000000+02 1.0171679-05 
9 1970 -1 12 9-999N935Z2401 1.0000000+02 1.8061980-06 
1 1100 -1 122 1.0000235+02 1.0000000+02 1.5391297-06 
iW 4110 “3 132 1.0000288+02 1.00000u0+02 7.6109598-06 
12 1120 -1 1u2 9. 9996854401 1.00000C0+02 1.2817393-06 
3 1150 -1 152 9.999T405+01 1.00000U0+02 4&.3133787-06 
Ly 11680 -\ 162 1.0000758+02 1.00000C0+02 2.2412374-05 
5% 1450 =.) 172 9.9989356+01 1.0000000+02 3.8367009-06 
16 1169 1 1769 1.0000774+00 1.0000000+00 6.4961967-02 
wv Vi70 1 1609 9.9991272-01 1.00600u00+00 4.9294689-02 
1A 1180 ) 1vo9 9.9976323-01 1.00000L0+0U 8.5637964-02 
19 Vivo 5 1s 1602 1702 1802 1902 9.99766886-91 1.0000000+00 8.54 73059-02 
SUBSTITUTIONS .-~ & SUBSTISUTION SETS, OF 11 ITEMS, EACH. 
Loc. 1ST SUBST. 2ND SUBST. BRD SUBST. 4TH SUBST. 
“193 -8,6990TH4+05 =1.06355291+07 3.B4O791IZ415 5.16044 55e14 
-“y 6.90000U0+0o 8.6060000*90 1.00v0000¢01 7.00000u0+0U 
i $.7653234902 3. 7653234402 3. 7655234402 3.7653238#02 
4% $.345795200) &.7129203¢01 &.9695KT71+0) &.0865409401 
20 ¥.0000000+03 9.000000U493 ¥.000000u+03 9.00000U0+05 
ot &. 10900G0+02 2. 1609000002 1. 1000000902 3. 1000000402 
vi 1.0000000+02 3.9009000+02 &, OULOO0U+02 2.00000U0+02 
2! %. 106090002 2-10000uUuU902 1. 1y000Gu+02 3. 10000u0*02 
151 6. 000000407 6-0009000+02 9.0000000+02 7.0000000+02 
22 1.1250Ju0¢02 6-692B912901 & 9282037401 8-5562077T+01 
23 2. 2600060902 1.75000CG+02 1.2500000+02 2.0000040+02 


Figure 6. 


Printed output showing final 
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prescription of zoom lens. 


Shows selected optical characteristics o: 


Positions. 
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A igwal Pe. Light resolut 


average is about 0.013 mm. 


of about 100 lines/mm is predicted for 


Zoom Element Focal 
Position Shift, 
mm 
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CONCLUSION 


This code has demonstrated its ability to analyze the 
performance of and/or to design simple or complex lens systems 
rapidly and precisely. Correlation between the lens perform- 
ance predicted by the program's multiple ray-tracing analysis 
and the measured performance of the constructed lens has been 
close. When there has been an apparent discrepancy between 
measured and predicted performance, the cause has always been 
a manufacturing error. It is therefore necessary to give the 
manufacturer construction tolerances of sufficient tightness 
to insure that design specifications will be met. The con- 
struction-tolerance specifications can be determined by the 
program's ability to evaluate the effect of small parameter 
changes and the displacement and/or tilt of the surfaces. In 
this connection, particularly for complex designs, one needs 
precise values for the refractive indexes of the glasses (melts) 
used for making the lens elements, since the manufacturing 
tolerance for glass stock is greater than the tolerance allow- 
able if predicted lens performance is to be obtained. 

In conclusion, the code gives to problems of lens design a 
general automated solution which makes full use of the digital 
computer's capabilities and abandons the classical approach 


requiring direction at every stage by the skilled specialist. 
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